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EFFECT OF MODIFIED THERMAL CONDUCTIVITY ON THE 
TEMPERATURE DISTRIBUTION IN THE PROTONOSPHERE 

by 

H. G. Mayr, Ernest G. Fontheim, and K. K. Mahajan 
ABSTRACT 

At typical protonospheric electron densities the electron mean free path is 
sufficiently long so that the coefficient of thermal conductivity is no longer 
given by Spitzer's expression. The effect on the temperature profile of using 
the corrected expression for conductivity is investigated. The corrected thermal 
conduction coefficient is density-dependent and has a more complicated tempera- 
ture dependence than the coefficient applicable to higher density plasmas. The 
results indicate that the effect is not negligible even under quiet conditions and 
at low latitudes. 
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EFFECT OF MODIFIED THERMAL CONDUCTIVITY ON THE 
TEMPERATURE DISTRIBUTION IN THE PROTONOSPHERE 


Recent daytime temperature measurements by two satellites at different 
points along the same magnetic field line inside the plasmasphere have allowed 
us to deduce the field aligned temperature gradients (see Fig. 1). It is, there- 
fore, of interest to investigate whether the presently used value of the heat 
conduction coefficient is consistent with these data. Recently, Mayr et al. (1972) 
have solved the simultaneous continuity, momentum, and energy equations for a 
multi-component ionosphere. Using their model we have investigated the effect 
of a modified coefficient of heat conduction on the temperature profile. It is 
well known that the thermal conductivity coefficient for plasmas obtained by 
Spitzer and Harm (1953) is no longer valid when the mean free path k approaches 
the characteristic length Z c given by (l/T dT/ds) -1 or even becomes larger 
than Z . Mayr and Volland (1968) have derived an improved expression for the 
thermal conductivity in the form of a power series in \/£ c . These authors define 
two mean free paths, k_ referring to particles approaching the point s = 0 from 
below and X + referring to those approaching from above, s = 0 is the point about 
which the expansion is taken. It is shown that to second order in A./'C the heat 
conduction flux can be expressed by the following relation: 


k _ . . dT 3k _ , x2 . 2 . d 2 T 

1 = c (k. + k ) n c (A/ - A“ ) n 

J * ' + ' e ds 32 + ' e ds 2 


( 1 ) 


K .. T S/2 *1. 
eff ds 


(2) 
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where c is the mean speed of the electrons, n g the electron density, k is Boltz- 
mann's constant, s the coordinate along the magnetic field line, and the effective 
conductivity k is defined by Equation (2). The conductivity is now a function 
of the electron density. 

If \ + 4 K , then \ > 0, because \ + refers to a higher altitude than 

does . Furthermore, the electron temperature has a maximum in the equa- 
torial plane (under symmetrical conditions), and therefore its second derivative 
is negative. As a result, the correction term in Equation (1) has the effect of 
reducing the thermal conductivity, as is to be expected. 

The mean free path has the general form 

T 2 

\ = K-L, (3) 

n 

where the constant K has been chosen such that * eff reduces to the conductivity 
of Spitzer and Harm if A. + = and the electron distribution is Maxwellian. 

In that case 

K ~ 3 x 10 4 cgs units. 

\ + and k_ in Equation (1) have been calculated from Equation (3) by using the 
values of T and n at distances k/2 above and below the point s = 0 respectively. 

In order to estimate the range of values of k „ along the field lines of 
interest here, the effective conductivity has been calculated for a number of 
different density profiles and a fixed temperature profile. The results are 
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shown in Figure 2. The curves are labelled by the electron density at 2500 km 
(equatorial plane). The values of k f{ in the shaded region correspond to 
realistic density profiles. The straight line represents the conductivity /< SH 
of Spitzer and Harm (1953) which is density-independent. Hence, it is seen that 
the correction to k is at most 15% in this particular case. 

SH 

Figure 3 shows two calculated temperature profiles along the 30° field 
line (equatorial radius = 2500 km). In addition the Explorer 22 and 32 tempera- 
ture measurements at 1000 km and 2500 km respectively along the same field 
line are shown. The rectangular boxes indicate midday ranges of altitude and 
temperature (see Fig. 1). These profiles have been obtained by solving the 
continuity, momentum, and energy equation’s so that each temperature profile is 
consistent with the electron density profile used in the energy equation. The 
lower boundary value of the electron temperature is T eo = T n = 1050°K. Curve 
(1) was calculated by using the conductivity k and a photoelectron flux at 
300 km of 2.5 x 10 9 cm -2 sec -1 'of 10 ev. The temperature gradient is clearly 
too small. 

If the modified heat flux given by Equations (1) and (2) is used, then one 
obtains curve (2). k cc was calculated with the temperature and density values 
corresponding to curve (1). The same photoelectron flux as in curve (1) was 
used to obtain the heating rates. The 0 + , H + , and electron density profiles 
corresponding to curve (2) are shown in Figure 4. These profiles were obtained 
simultaneously with the temperature profile and are based on an upward H + 
flux in the protonosphere of 2.7 x 10 7 cm" 2 sec" 1 , typical of daytime conditions. 
This flux value was chosen in order to obtain electron density values in agree- 
ment with the measured ones. 
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The effect of the diminished conductivity is to increase the temperature 
gradient, as is to be expected. However, the calculated gradient is still smaller 
than the observed one. it should also be pointed out that this calculation was 
carried out for a relatively low latitude. Since the corrected thermal conductivity 
is density-dependent and the density decreases toward higher L-values, it is 
clear that at higher L-values the effect on the temperature profile due to 
the corrected conductivity coefficient will be even more pronounced. 
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Figure 1. Combined T g measurements for Explorers 32 and 22 along a single 
field tube L = 1.4. The difference in temperature at these two 
altitudes reflects the temperature gradient along this tube (from 
Brace, 1970). 

Figure 2. Calculated effective heat conductivity for different density profiles. 

Figure 3. Calculated temperature profiles for the cases of the conductivities 
/< gH (curve 1) and * eff (curve 2). 

Figure 4. 0 + , H + , and electron density profiles corresponding to the tempera- 

ture profile of curve 2 in Figure 3. 
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Figure 1. Combined T e measurements for Explorers 32 and 22 along a single field tube L =1.4. 
The difference in temperature at these two altitudes reflect the temperature gradient along this 
tube (from Brace, 1970). 
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Figure 4. 0 + , H + , and electron density profiles corresponding to the temperature 

profile of curve 2 in Figure 3. 
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